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Abstract. The paper studies a non-linear transformation between Brown- 
I ian martingales, which is given by the inverse of the pricing operator in the 

mathematical finance terminology. Subsequently, the solvability of systems 

C^l of equations corresponding to such transformations is investigated. The lat- 

ter give rise to novel monotone pathwise couplings of an arbitrary number 
of certain diffusion processes with varying diffusion coefficients. In the case 
that there is an uncountable number of these diffusion processes and that 
the index set is an interval such couplings can be viewed as models for the 
growth of one-dimensional random surfaces. With this motivation in mind, 
we derive the appropriate stochastic partial differential equations for the 

i— i growth of such surfaces. 

Oh 

> 

1. Introduction 
The starting point of the paper is the equation 

> (1.1) M(t)=E[v(X(T))\F x (t)], te[0,T\, 

CN| where T is a positive real number, M(t), t e [0, T) and X(t), t e [0, T] are 

stochastic processes with continuous real-valued paths, (J rX (t)) te [ 0i T] is the fil- 
tration generated by X and v : R — > R is a (deterministic) function. This type 
^ of equations is fundamental to the mathematical theory of asset pricing, where 

y— i the process X stands for the price process of an asset in a financial market, 

v(X(T)) for the payoff of a contract on this asset and M for the price process 
of this contract. In this context, one usually assumes that the process X is of 
a particular form (for example, a logarithmic Brownian motion in the classical 
Black-Scholes model as in [3], [TSj) and tries to determine exactly or to obtain 
estimates on the process M. It is standard to model the asset price process 
X by a continuous diffusion process and one assumes that X is a martingale, 
which implies the absence of arbitrage in the market given by X and M by the 
Fundamental Theorem of Asset Pricing (see [9], pU], [6], [7] and the references 
therein). For the rest of the paper, we will call continuous diffusion processes, 
which are martingales, Brownian martingales. 
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The starting point of our analysis is the following question asked by David 
Aldous ( p] ) : Suppose that M is a given Brownian martingale and that the func- 
tion v is given. Can one find a process X in the class of Brownian martingales 



such that (1.1) holds in the sense of equality of probability laws? In other words, 
can one find a Brownian martingale X such that the process K[v(X (T))\J rX (t)], 
t G [0, T] coincides in law with a given Brownian martingale Ml The first main 
result of the paper is that, under certain regularity assumptions on M and v, 



the question can be answered in the affirmative. The operation of solving ( 1.1[ ) 



for X can be viewed as the inverse of the pricing operator, which maps (X, v) 



to M via (1.1). We remark at this point that, although we allow M and X to 
take negative values throughout, one can enforce the nonnegativity of M and 
X suggested by the financial interpretation by letting M be a Brownian mar- 
tingale taking only nonnegative values and choosing v such that the preimage 
of the state space of M under v is a subset of [0, oo). Similarly, one can ensure 
that M and X take values in the desired bounded subintervals of [0, oo). 

In the second part of the paper, we study systems of equations 

(1.2) M A (-) =E[v x {X{T))\F x {-)], AG A 

for arbitrary index sets A. We provide a necessary and sufficient condition on 
the family M A , A G A of Brownian martingales under which it is possible to find 
a Brownian martingale X and functions v x : R — > R, A G A of certain regularity 



satisfying (1.2). Moreover, if A is endowed with a partial order -<, we give a 



condition for the system ( 1.2 ) to be solvable in the specified sense with functions 
v x : R — y R, A G A, which satisfy v Al < v x ' 2 whenever Ai -< A2. Clearly, the 
result are monotone pathwise couplings of the Brownian martingales M A , A G A. 
If A C R is an interval endowed with the usual total order, then by viewing A 
as the time variable and t as the space variable, one may regard the process 
M x (t), (A, t) G A x [0, T] as a model for the growth of a one- dimensional random 
surface over time. We show that the growth of the surface can be described by 
a stochastic partial differential equation (SPDE). 

To state our first result rigorously, we introduce some notation. In the context 



of equation (1.1), we denote the state space of M by S C R, which we assume 



to be a (bounded or unbounded) open interval, and write C t = \ a(t, x) -4^ for 
the time-dependent generator of the process M. Moreover, we let Q be the 
space of continuously differentiable functions v : M. — > R, whose derivatives v' 
are bounded between two positive constants. 

Theorem 1. Suppose that S is contained in the range of v and that the tran- 
sition densities p M (t,x;T,y), t G [0,T), x,y G S of M exist and possess con- 
tinuous partial derivatives pf 1 , p^f , p^f x and p^ xx . Moreover, assume that there 
are constants C\,C<i > such that 



1.3) / \pM(t,x;T,y)\(l + \y\)dy <C x e c ^ 

Js 
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for all (t, x) G [0, T) x S and that for any given t G [0, T) and x £ S there exists 
an open neighborhood U of (t, x) and a function F : S — > [0, oo) such that 

(bf I + bf I + bf I + ]rfL\)(t,2;T,y) < F(y), (t,x,y) e U x S, 



[ F(y)(l + \y\)dy<oo. 
Js 



Finally, suppose that a takes only positive values, its partial derivative a x is 
continuous and the estimates 

(1.4) \a(t,x)\ < C 3 , \a x (t,x)\ < C 3 (l + |x|), (t,x) G [0,T] x 5 

hold for some constant C3 > 0. Then, for any v G Q, there exists a unique 
Brownian martingale X satisfying the equation (1.1) in law and it solves the 
stochastic differential equation (SDE) 

(1.5) dX(t) = h x (t, /i ( - 1) (t,X(t)))y / a(t,/i(- 1 )(t,X(t)))d J B(t), t G [0,T], 
where h solves the backward Cauchy problem 

(1.6) ^ + ^^ = 0, fc(r, •) = ^ 

in i/ie classical sense, the superscript (—1) denotes the spatial inverse and B is 
a one- dimensional standard Brownian motion. 



Our main result on the solvability of the system (1.2) requires the following 
definition. 

Definition 1. Let M A , A G A be a family of Brownian martingales with time- 
dependent generators \ aX (t, x ) 2 -^2 and initial values mj, respectively, and fix 
the notations £ A (t,x) = A G A and A = (S A )(- 1 ), A G A. We 

call the family M A , A G A consistent if the following conditions hold. For any 
Ai 7^ A2 in A, there is a function & Al > A2 : [0, T] — >• R, which solves the ordinary 
differential equation (ODE) 

= — a ^{t,Q Xl (t^ M + 6 Al ' A2 )) + S Al (t,0 Al (t,S A2 +6 Al ' A2 )) 

(1.7) 2 

+-a A2 -£ A2 . 

for all values of x such that the pointwise limit 

(1.8) r Al ' Aa := lim© Al (t,S A2 + b XlM ) 

exists, belongs to Q, and one has 

(1.9) m Aa = E' m o 1 [(^"^-^(M^T))]. 

Hereby, the argument of 6 Al,A2 , 6 Al,A2 is £, and the arguments of S A2 , er A2 , S A2 
are t, x. 

It turns out that this notion of consistency is necessary and sufficient for the 
solvability of the system fll.2|). For a discussion of the conditions in Definition 
[TJ please see Remarks [T] and [2] following the statement of Theorem [2] 
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Theorem 2. Suppose that M x , A G A is a family of Brownian martingales, 
which all satisfy the conditions of Theorem [TJ In addition, suppose that the 
diffusion coefficients o~ x , A G A of M x , A G A are continuously differentiate 



and bounded away from 0. Then, the system (1.2) is solvable by a Brownian 
martingale X and a family v x , A G A of functions in Q if and only if the 
family M x , A G A is consistent. Moreover, if this is the case, then there are 
uncountably infinitely many such solutions and the following are true: 

(a) The processes M x , A G A can be defined on the same probability space to 
form a weak solution of the degenerate system of SDEs 

(1.10) dM x (t) = a x (t, M x {t)) dB(t), M A (0) = m A 
on [0, T] . 

(b) If, in addition, there is a partial order -< on the set A and the functions 
pAi,A2 i n Definition^ are such that Xi ~< X 2 implies x > r Al ' A2 (x) for all 
x in the state space of M X2 , then the processes M x , A G A can be defined 
on the same probability space in such a way that, for all Ai -< A2 in A, the 
inequality M Xl (t) < M X2 (t) holds for all t G [0,T] with probability 1. 

(c) If, in the situation of part (b), the set A C R is an interval with the usual 
total order, the function A 1— > m A is continuously differentiate and the 
function o~(X, t, x) := o~ x (t, x) has continuous and bounded partial derivatives 
o~\, a x\, o~xx, o~x and o xx , then there is a version of the growth process 
H(X,t) := M x (t), (X,t) G A x [0,T] ; which is continuously differentiate in 
X, and such that the corresponding derivative if a solves 

(1.11) H x = (a x (X, t, H) + ^.(A, t, H)H X ) £, 

(1.12) fT A (0,A) = ^. 

Hereby, ^ is a distribution-valued Gaussian field with correlation function 

(1.13) E[f (ti, Ai)£(t 2 , A 2 )] = 6{t 2 - h), Ai, A 2 G A, t h t 2 G [0, T]. 



Remark 1. Consider first the situation |A| = 2 and let M be a Brownian mar- 
tingale with diffusion coefficient a and initial value mo, which satisfies the 
conditions of Theorem [TJ Then, a family of Brownian martingales with diffu- 
sion coefficients a and initial values mo, which are consistent with M, can be 
constructed as follows. First, let be a solution of the backward heat equation 

Q t + \ e** z^Ht, x + d)) + E t (t, ^-'\t, x + D))j e x = 

with some terminal condition 0(T, •) G Q, where D is an arbitrary real constant. 
Next, set 

(1.14) S = 6 ( " 1) , a = l/t x and m = E m ° [S ( - 1} (T, S(T, M(T)) + D)}. 
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Figure 1 . The figure demonstrates the growth of a random sur- 
face as in Theorem [2] (c). The family of consistent Brownian 
martingales is the one of Example [T] at the end of subsection 3J2 
and the five curves correspond to sample paths of the growth 
process if for A = 1, 2, 3, 4, 5. 



Then, the Brownian martingale M corresponding to the diffusion coefficient a 
and initial value mo is consistent with M. Indeed, by differentiating the equa- 
tion 6(t, x)) = x once with respect to t and once and twice with respect 
to x, one can express the partial derivatives of in terms of the partial deriva- 
tives of S. Plugging the resulting formulas into the backward heat equation and 
using a = one then easily verifies that the constant function b(t) = D 

solves the ODE in Definition [T] 

In the case that a = 1, that is, when M is a standard Brownian motion, the 
backward heat equation simplifies to 

(1.15) o = e t + ^e M . 

Next, take 0(T, •) to be an arbitrary function in Q. Then, 0(T, •) can be 
viewed as the cumulative distribution function of an infinite positive measure 
fi on M.. Clearly, for any t G [0,T), the function 0(t, •) is given by a cumulative 
distribution function of the measure fi * <pT-t, where <pT-t is the normal density 
with mean and variance T-t. Consequently, E(t, •) = e (-1) (t,-) is given 
by a quantile function qr-t of // * <pT-t- Thus, we conclude that the Brownian 
martingale solving 

(1.16) dM(t) = — I— dB(t) 

q' T _ t {M{t)) 
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is consistent with a Brownian motion started in m . We note hereby that, since 
we can choose the constant D above in an arbitrary manner, we can take mo 
to be an arbitrary real number. 

Finally, in the case that A has more than two elements, it suffices to observe 
that Theorems [T] and [2] imply that the definition of consistency for pairs of 
Brownian martingales induces a transitive relation on the space of Brownian 
martingales satisfying the conditions of Theorem [1} Indeed, if M 1 , M 2 , M 3 are 
Brownian martingales such that M 1 and M 2 are consistent and M 2 and M 3 
are consistent, then by Theorem |2j there are functions v 1 ,v 2 ,w 2 ,w 3 G Q and 
Brownian martingales X, Y such that 

M\.) = E[v\X(T))\J^ x (.)], M 2 (.) = E[v 2 (X(T))\F x (.)}, 
M 2 (-) = E[w 2 (Y(T))\F Y (-)}, M 3 (-) = E[w 3 (Y(T))\F Y (-)}. 

Now, the uniqueness statement of Theorem [T] implies that the processes X and 
(t> 2 )(~ 1 )(?/; 2 (y(-))) coincide in law. Therefore, the process 

E^^^V^))))^*)], te[0,T\ 

has the law of M 1 , so that M 1 and M 3 are consistent by Theorem [2] This shows 
the transitivity of the consistency property for pairs of Brownian martingales 
and, thus, that the family of processes M constructed above is consistent, pro- 
vided that they all satisfy the conditions of Theorem [T] 



Remark 2. The following is the simplest example of a situation, in which the 
conditions of Definition [T] hold. Suppose that for every A G A the diffusion 
coefficient a x of M x does not depend on t, and that for all Ai 7^ A 2 in A there 
are constants A XlM > 0, B XlM G R such that 

a Xl (x) = -j^a X2 (A Xl ' X2 x + B XlM ) 

and tjiq 2 = A XlM m Xl + B XlM . Then, the family M A , A G A is consistent. 
Indeed, choosing 6 = and r Al ' A2 (x) = x ~ A f - 2 , one can check that the con- 
ditions of Definition 1 are satisfied by using the identities S A ^. = — ^\^ 2 an d 
y\-2 = a? 

^XX (cr A 2)2- 



Remark 3. Considering the process H in Theorem [2] (c) as given, one can view 
the SPDE (1.11) for the partial derivative H\ as a degenerate version of the 
SPDEs studied in [3], [5], with random instead of deterministic coefficients. 
Indeed, setting A = in equation (1.1) of |4J and choosing the coefficient g\ 
there to be the (random) diffusion coefficient in our equation (1.11) and the 
coefficient / there to be the corresponding (random) Ito correction term, one 
recovers our SPDE (1.11). For other SPDEs with noise, which is white in time 
and correlated in space, we refer the reader to [15], [IE] and the references 
therein. 
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The rest of the paper is structured as follows. Section 2 is devoted to the 
proof of Theorem [TJ As it turns out, equation (1.1 ) can be solved by solving a 



nonlinear parabolic partial differential equation (PDE). This PDE is similar to 
the nonlinear PDEs appearing in [19] and can be reduced to a linear PDE by a 
suitable transformation. This is the content of subsection 2.1. In subsection 2.2, 
we solve the linear PDE and complete the proof of Theorem [TJ In subsection 



2.3, we analyze the solvability of equation (1.1) in the degenerate case that 
v — l[c,oo) for some cel. In mathematical finance terms, the equation (1.1) 
for this choice relates the price process M of a digital option to its payoff 
l[ coo )(X(T)). Alternatively, one can view the process M as describing the 
evolution of the conditional probabilities of the event {X(T) > c} based on 
the information about the process X available so far. Such processes can be 
observed in practice as quotes in prediction markets (see [2] for more details). 
In contrast to the findings of Theorem [TJ we find an uncountably infinite family 



of Brownian martingales solving the equation (1.1) in this case. 

In section 3, we give the proof of Theorem 
first prove the relevant statements of Theorem 



For the sake of clarity, we 
in the case that |A| =2. To 
this end, we need to solve a system of two nonlinear PDEs, each of the same 
type as in subsection 2.1. This is the content of subsection 3.1. In subsection 
3.2, we show how the arguments extend to the case of general index sets A. 
Finally, in section 4 we give more examples of situations, in which Theorems 
[TJ and [2] apply, by treating Brownian martingales consistent with the Kimura 
martingale. 

2. Proof of Theorem [TJ 
2.1. Reduction to a linear PDE. In this subsection, we show that the prob- 



lem of solving (1.1) naturally reduces to a backward Cauchy problem for a 



linear PDE. We start with a proposition. 

Proposition 3. Let the functions a and v be as in Theorem^ and suppose that 



X is a Brownian martingale solving (jl.l|) with this v and M being a Brownian 

dx 2 



martingale with time- dependent generator \a{t } x) 4^. Moreover, assume that 



there exists a classical solution u of 

1 aft, k) 

(2.1) u t + -^^u xx = 

such that u(T, •) = v, the partial derivative u xxx exists and is continuous, and 

(2.2) \u x \ > Cie _C2 '"' for some constants ci, C2 > 0. 

Then, the time- dependent generator of X is given by | a ^'"|^^ and it holds 

(2.3) E[v(X(T))\F x (t)} = u(t,X(t)), t e [0,T]. 
The proof of the proposition relies on the following lemma. 

Lemma 4. Let the functions a and v be as in Theorem [TJ and suppose that 



u is a classical solution of the PDE (2.1) satisfying u(T,-) = v and (2.2), 
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whose partial derivative u xxx exists and is continuous. Then, the spatial inverse 
h := w^ 1 -* is well-defined, solves the linear backward Cauchy problem 



(2.4) 



ht + - a(t, x) h a 



0, h(T, 



in the classical sense and its partial derivative h x is bounded. 



Proof. From the PDE (2.1) it is clear that either u x > everywhere, or u x < 
everywhere. However, since v is strictly increasing, it holds u x (T, •) = v' > 
and, hence, the inequality u x > must be true everywhere. Thus, for every 
t G [0,T], the function u(t,-) is strictly increasing and therefore the spatial 



inverse h 



u 



(-1) 5 



is well-defined. Next, we consider the equation u(t, h(t,x)) 



x. Differentiating this equation with respect to t once and with respect to x 
once and twice, we obtain the following set of equations: 

(2.5) ut(t, h(t, x)) + u x (t, h(t, x))h t (t, x) = 0, 

(2.6) u x (t, h(t, x))h x (t, x) = 1, 

(2.7) u xx (t, h(t, x))h x (t, x) 2 + u x (t, h(t, x))h xx (t, x) = 0. 

Hereby, the partial derivatives ht, h x and h xx exist by the Implicit Function 
Theorem. Plugging the equations (2.5), (2.6), ( 2.7| ) into (2.1), we obtain 



{21 



h t (t,x)u x (t, h(t,x)) + 



a(t,x) -u x {t,h{t,x))h xx {t,x) 



2(l/h x (t,x)y 



h x (t,x) 



0. 



Simplifying, we conclude that h is a classical solution of the problem (2.4) 



as desired. Moreover, differentiating both sides of the equations in (2.4) with 
respect to x, we deduce that h x is a classical solution of the problem 



(2.9) (h x ) t + -a x {t,x){h 



rlx + -a{t,x) (h a 



0, h x (T, 



Hereby, the existence and continuity of the partial derivative h xxx follows from 
the existence and continuity of u xxx and the Implicit Function Theorem. More- 
over, the existence and continuity of the partial derivative h tx is a consequence 

In addition, by as- 
is bounded. 



of (2.4) and the existence and continuity of a x and h a 



sumption, the function v' is bounded away from 0, so that h x (T, 
Moreover, in view of (2.6), the assumption (2.2) implies 



(2.10) 



\hn 



1 

< — e 

Cl 



c 2 M 2 



Applying the Maximum Principle for linear parabolic PDEs in the form of 
Theorem 9 in chapter 2 of [8J , we conclude that h x must be bounded everywhere. 
This finishes the proof of the lemma. □ 



Proof of Proposition [3[ B y Lemma |4| the spatial inverse h = ^ is well- 
defined, solves the PDE ( |2.4| ) in the classical sense and has a bounded partial 
derivative h x . Defining the process N(t) := E[v (X(T))\F X {t)], t E [0,T], we 
obtain by Ito's formula: 

(2.11) dh(t,N(t)) = h x (t,N(t))dN(t), te(0,T). 
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Hereby, the drift terms disappear, since N has the time-dependent generator 
| a(t, x) 4-$ by assumption and h is a classical solution of the equation (2.4) 



Moreover, since h x is bounded and N is a martingale, the process h(t,N(t)), 
t G [0,T] is also a martingale. Noting that h(T,N(T)) = u^ l \T,v(X{T))) = 
X(T) and recalling that X is a martingale, we conclude that h(t,N(t)) = X(t) 
must hold for all t G [0, T] with probability 1. Clearly, this implies that N(t) = 



u(t,X(t)), t G [0, T] with probability 1. Hence, from equation (2.11), we see 
that the time-dependent generator of X must be given by 

1 , . . . , 9 . . d 2 1 a(t,u(t,x)) d 2 
(2.12) -h x {t,u{t,x)) 2 a{t,u{t,x)) 



2 ' ' ' ' dx 2 2 u x (t,x) 2 dx 2 

Finally, by the definition of N and the previous considerations, it holds 
(2.13) E[v(X(T))\F x (t)\ = N(t) = u(t,X(t)), t G [0,T]. 

This finishes the proof of the proposition. □ 



The proposition shows that, if a classical solution u to the PDE (2.1) with 
terminal condition v exists and possesses a continuous partial derivative u xxx , 



then a solution X of the equation ( 1.1 ) must have the time-dependent generator 



2 ^ujt^ ^ n ^ e nex ^ subsection, we show that, under the assumptions of 



Theorem [TJ the desired classical solution of the nonlinear PDE (2.1) exists by 



solving the linear PDE (2.4), and that the resulting process X is indeed a 
Brownian martingale. 

2.2. Solution of the linear PDE. As we have seen in the previous subsec- 



tion, the crucial step in the solution of the equation (1.1) consists of finding a 



sufficiently regular solution of the problem (2.4). In the following lemma we 
show that the latter exists under the assumptions of Theorem [TJ 

Lemma 5. Let the assumptions of Theorem^ hold. Then, the problem (2.4) 
has a classical solution h, which satisfies 

1 I 1 2 

(2.14) \h x \ < -e C2m for some constants ci, c-i > 

Cl 

and whose partial derivative h xxx exists and is continuous. 

Proof. Recalling the notation p M (t,x;T,y), t G [0,T), x,y G S for the transi- 
tion densities of the process M, we define the function h by 

(2.15) h(t,x) = [ p M (t,x;T,y)v(- 1 \y)dy. 



s 



The assumptions on p allow us to interchange the order of differentiation and 
integration to obtain 



(2.16) h t (t,x) = / pf I (t,x-T,y)v ( - 1 \y)dy, 

(2.17) h xx (t,x) = [ p x 4 x (t,x;T,y)v(- 1 \y)dy, 



s 
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and to conclude that the functions h t and h xx are continuous by using the 
Dominated Convergence Theorem. It follows that 

tk{t,x) + ^a(t,x) h xx (t,x) 

[pf(t, x; T,y)+ l - a(t, x) p%(t, x; T, y)] v^(y) dy = 



is a consequence of the Kolmogorov backward equation 
(2.18) pf(t, x- T,y)+ l - a(t, x) p™(t, x; T, y) 



0. 



Moreover, by the definition of h, we have h(T, •) = ?/ 1 \ In addition, the growth 



estimate (2.14) follows from another exchange of the order of differentiation and 
integration 



(2.19) 



h x (t,x) = / pM(t,x-T,y)v ( - 1 \y)dy, 



G Q and (1.3). Finally, the partial derivative h xxx exists and is continuous 



due to yet another exchange of the order of differentiation and integration 

(2.20) K xx (t,x) = [ p x l x (t,x;T,y)v ( -- 1 \y)dy 

Js 

and the continuity of the right-hand side in the latter equation, justified by the 
assumptions on p M and the Dominated Convergence Theorem. □ 

We are now ready to give the proof of Theorem [TJ 



Proof of Theorem [TJ We first prove that a solution X to ( 1.1 ) of the described 
form exists. To this end, we let h be a classical solution of the problem (2.4) 
as in Lemma|5]and set X(t) = h(t, M(t)), t G [0, T]. Applying Ito's fomula, we 
see that 

(2.21) dX(t) = h x (t,M(t))dM(t), te[0,T]. 

Hereby, the drift terms disappear, because M has the time-dependent generator 



\ a(t,x) 4^2 and h is a classical solution of the PDE in (2.4). Moreover, noting 
that h x is a classical solution of the problem (2.9) and applying the Maximum 



Principle for linear parabolic equations in the form of Theorem 9 in chapter 2 of 



[H] (note the growth estimate (2.14)), we conclude that h x is bounded between 
two positive constants. In particular, the spatial inverse u := h^~^ is well- 
defined and M{t) = u(t,X(t)), t e [0,T]. Putting these observations together 



with (2.21 ), we conclude that X is a Brownian martingale with time-dependent 



generator \ h x (t, u(t, x)) 2 a(t, u(t, x)) 

Next, we note that the identities X{t) = h(t,M(t)), M(t) = u(t,X(t)), 
t G [0,T] imply that the filtration generated by M is the same as the filtration 
generated by X. Moreover, we have 

(2.22) v(X(T)) = v(v { - 1] (M(T))) = M(T). 
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Hence, since M is a martingale in the filtration generated by X, it holds 

(2.23) M(t) =E[M(T)\F x (t)} = E[v(X(T))\F x (t)\, t E [0, T] 
as desired. 

We now turn to the proof of uniqueness. Differentiating both sides of the 
equation h(t, u(t, x)) = x once with respect to t and twice with respect to x, 
we get 

(2.24) h t (t, u(t, x)) + h x (t, u(t, x))u t (t, x) = 0, 

(2.25) h xx (t, u(t, x))u x (t, xf + h x (t, u(t, x))u xx (t, x) = 0. 



Hence, 



1 a(t, u) h t (t, u) + la(t, u)h xx (t, u) _ 



In addition, the existence and continuity of the partial derivative h xxx and the 
Implicit Function Theorem imply that u xxx exists and is continuous. More- 



over, the estimate \2.2\ is a direct consequence of the growth estimate (2.14) 




Applying Proposition pi we conclude that a Brownian martingale X solving 



(1.1) must have the generator | °„ '^fi^ This and the boundedness of the 



function u x show that the process u(t, X(t)), t G [0,T] is a martingale in the 



filtration generated by X by an application of Ito's formula. Therefore, by (1.1 ) 
(2.26) u(t,X(t)) = E[v(X(T))\F x (t)} = N(t), t 6 [0,T] 



for a process iV of the same law as M. From (2.26) we see that the law of the 



process X(t) = h(t,N(t)), t 6 [0, T] is uniquely determined. □ 

2.3. A degenerate case. We turn now to the degenerate case that the func- 
tion v in equation ( 1.1[) is given by the indicator function l[ c ,oo) f° r some c G 



To make sure that (1.1 ) is well-posed, we need to assume that the state space of 



M is the interval [0, 1] and that M(T) G {0, 1} holds with probability 1. As be- 
fore, we write | a(t, x) 4^ for the time-dependent generator of M. In contrast 
to the non-degenerate case of Theorem [TJ there exist uncountably infinitely 



many Brownian martingales solving (1.1) here 



Proposition 6. In the case that v = l[ Cj0 o) for some c EM. in equation ( 1.1 ) and 



that M is a Brownian martingale as just described, there exists an uncountably 



infinite family of solutions X of equation (1.1), each of which satisfies an SDE 
of the form 

(2.27) dX{t) = (ci - Co) ( X W " c o)/(ci - c )) dB{t), t E [0,T], 
for some constants Cq < c < C\, where B is a standard Brownian motion. 

Proof. As in the non-degenerate case, we consider the linear PDE 

(2.28) h t + - a{t, x) h xx = 0, 
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which formally yields 

h(t,x) = E[h(T,M(T))\M(t) = x] = xh(T, 1) + (1 - x) h{T, 0). 

This motivates letting c\ := h(T, 1) > c, Co := /i(T, 0) < c and setting = 
c\M(t) + co(l — M(t)), t G [0,T]. Indeed, then X is a Brownian martingale. 
Moreover, 

E[l{x ( T)> c} |^ X (t)] = E[M(T)\F x {t)} = M(t), 
since the filtration generated by X is the same as the filtration generated by 
M, and M is a martingale in this filtration. It follows that X is a solution of 



the equation (1.1). Finally, by Theorem 4.2 in chapter 3 of [12] . the process M 
obeys the SDE 



(2.29) dM(t) = y/a(t, M(t)) dB(t), t G [0, T] 

for a standard Brownian motion B. This shows that the process X(t) = 
cxM(t) + Co(l — M(t)), t G [0, T] satisfies the SDE in the statement of the propo- 
sition. Now, it remains to note that the choice of the constants Cq < c < c\ was 
arbitrary and the proof is finished. □ 

3. Proof of Theorem |2] 

3.1. Systems of two equations. For the sake of clarity of exposition and sim- 
pler notation, we first give the proof of the relevant statements of Theorem [2] in 
the case that |A| = 2. To fix notations, we let M, M be two Brownian martin- 
gales with time-dependent generators C t = \ a(t, m) 2 -£^2, A = ^cr(^; m ) 2 ^2 
and initial values m , m , respectively. Assuming that M and M satisfy the con- 
ditions of Theorem [TJ we seek Brownian martingales X and functions v, v G G 
such that the system of equations 

(3.1) M(.)±E[v(X(T))\F x (.)}, 

(3.2) M(-) =E[d{X{T))\F x {-)] 

is satisfied. In the case that |A| = 2, the definition of consistency (Definition [l] 
in the introduction) simplifies to the following definition. 

Definition 2. Letting £(t, x) = J^ o and S(t,x) = JT q a7ty)> we sa y ^ na ^ 
the Brownian martingales M and M are consistent if there is a function b : 
[0, T] E, which solves the ODE 

(3.3) b = - l -a x {t, S ( " 1} (t, t + b)) + E t (t, ^- l \t, E + b)) + h x - % 
for all values of x and such that the pointwise limit 

(3.4) r(m) :=limE (-1) (t,I! + 6) 

exists, belongs to Q, and it holds 

(3.5) m = E mo [T(- 1 \M(T))}. 

Hereby, the argument of b, b is t, and the arguments of E, <j x , E t are t, x. 
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We show now that the system (3.1 ), (3.2) is solvable if and only if the Brow- 



nian martingales M and M are consistent. 

Proposition 7. Suppose that the Brownian martingales M, M satisfy the as- 
sumptions of Theorem [T] and that their diffusion coefficients a, a are contin- 
uously differentiable and bounded away from 0. Then, there is a Brownian 
martingale X and functions v,v £ Q, which solve the system ( 3.1[ ), (3.2), if 
and only if M and M are consistent in the sense of Definition |2 Moreover, if 
this is the case, then there are uncountably infinitely many such solutions. 



Proof. We assume first that M and M are consistent in the sense of Definition 
[2j In order to solve the system (3.1), (3.2), it suffices to find a Brownian mar- 
tingale X with a time- dependent generator ^a x (t, x) 
v,v £ Q such that the system of PDEs 



4-2, as well as functions 



(3.6) 


1 x 

u t + -a u xx = 0, 


(3.7) 


1 x 

u t + -a u xx = 0, 


(3.8) 


u x \fa~x = a(t, u), 


(3.9) 


u x \fa* = a(t, u), 


(3.10) 


u ( T i •) = v , 


(3.11) 


u(T, -) = v 



has a classical solution, E[t>(X(T))] = m and E[t>(X(T))] = m . Indeed, then 
the processes u(t,X(t)), t £ [0, T] and u(t,X(t)), t £ [0,T] satisfy 



(3.12) 
(3.13) 



du(t,X(t)) =a(t,u{t,X(t)))dB(t), te [0,7], 
du{t, X(t)) = a(t, u(t,X{t))) dB(t), t £ [0, T] 



for a standard Brownian motion B (apply Ito's formula and ( 3.6[ ), (3.8), (3.7), 
(3.9)) and, thus, are Brownian martingales with time-dependent generators 
C t = | a(t, m) 2 £^ , l t = \ &(t, m 
respectively. Since 



^ 2 ^2 and initial values u(0, X(0)), u(0,X(0)), 



E[v(X(7))\F x (.)] 
E[v(X(T))\F X {')) 



E[v(X(T))} 
E[v(X(T))} 



m , 
m , 



this gives the desired solution of the system (3.1), (3.2). Note hereby that by 
our assumptions on a and a, the solutions to (3.12) and (3.13) are pathwise 
unique and, therefore, the martingale problems for C t and C t are well-posed 
due to the results of Yamada and Watanabe (see e.g. Proposition 3.20 in [12J). 



To solve the system (3.6 )-( 3. 11 ), we fix an arbitrary function v £ Q and 
proceed as in the proof of Theorem [T] to find a Brownian martingale X and a 
function u such that the equations (|3.6[), (3.8), (3.10) and mo = E[v(X(T))] 
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are satisfied. Next, to make sure that equation (3.9) holds, it suffices to choose 
u in such a way that 

(3-14) = 

<J{t,U) <J{t,U) 

However, recalling the definitions of the functions S and £ in Definition [2] and 



integrating both sides of the equation (3.14) in x, we see that the latter is 
equivalent to the equation 

(3.15) £(*,«) -£(t,u) = b(t), te[0,T] 

for some function b : [0,T] — > R. Moreover, since the functions a and a take 
only positive values by assumption, the functions E(t, •) and £(£, •) are strictly 



increasing for all t £ [0, T\. It is now easy to check that equation (3.15) will 
hold if and only if we set u = hS~ 1 ^ with a function h satisfying 

(3.16) h(t,x) =M ( " 1) (t,S ( - 1) (t,S(t,x) + &(*))) =: u (-1) (t,r(t,x,6(t))). 

Hereby, we recall from the proof of Theorem [l] that the functions u(t, •), t £ 
[0, T] are strictly increasing, so that the spatial inverses ^'^(t, ■), t £ [0, T] are 
well-defined and strictly increasing as well. We note at this point that the just 
introduced function r is determined by the parameters of the problem. On the 
other hand, we are free to pick a function b of our choice and choose it as the 
function b in Definition [2] 



Next, we make sure that equation ( 3. 11[ ) holds. To this end, we note that by 



the definition of the function T in Definition |2| and (3.16) we have 
(3.17) h(T,.)=u^ 1 \T,T(.))=v^ 1 \T(.)). 
Thus, to ensure (3.11) it is enough to set v = I^ -1 ) o v. The resulting function 



v belongs to Q by our assumptions on T and v. Moreover, one has 
(3.18) E[v{X{T))\ = E[r(- 1 )(t;(X(T)))] = E mo [T^ (M(T))] = m , 



where the second identity is a consequence of the fact that X, v solve (3.1) 
by construction and the third identity is a consequence of M and M being 
consistent (see Definition [2]). 

It remains to check that with the choice of h (and, hence, also of u) above, 



the equation (3.7) is satisfied. We claim that it suffices to show that h is a 



classical solution of the PDE 

(3.19) h t +-a(t,x) 2 h xx = 0. 

Indeed, differentiating the equation h(t,u(t,x)) = x once with respect to t and 
once and twice with respect to x, one can express the partial derivatives of u 
in terms of the partial derivatives of h and compute 

~ 1 x~ — ' 1 1 dr(t,u) 2 _ h t (t,u) 1 a(t,u) 2 h xx (t,u) 

2 2 ul hx (t,u) 2 hx ( t ,u) 
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where we have used (3.9) in the first identity. To check that h is a classical 
solution of (3.19), we set h = u^ 1 ^ and observe that, in view of (3.16), equation 

d 2 



(3.19) can be rewritten as 
d 



1 



(3.20) ^h(t, r(t, x, b(t))) + | a(t, x) 2 — 



h(t,r{t,x,b(t))) = 0. 



Simplifying the left-hand side, we see that this is equivalent to the following 
ODE for b: 



(3.21; 



b(t) 



la 



(t, x) 2 (h xx r 2 x + h x r xx ) + h t + h x r t 

n h x 



Hereby, the arguments of h and its partial derivatives are t, r(t, x, b(t)), and the 
arguments of r and its partial derivatives are t, x, b(t). The key observation 
is now that, due to the definitions of the functions r, £ and £ (see (3.16) and 
Definition [2]), we have 

a(t,x) 2 r x {t,x,b) 2 = a(t, x) 2 £i _1) (t, £(t, x) + b) 2 t x (t, x) 2 
= E^t, t(t, x) + b) 2 = a(t, r(t, x, b)) 2 . 

Plugging this into (3.21 ) and recalling from the proof of Theorem |T| th at h is a 
classical solution of the problem (2.4), we can simplify the ODE (3.21) to 



i % 



(3.22) 



b(t) 



a{t, x) 2 r xx (t, x, b(t)) + r t (t, x, b(t)) 



r b (t,x,b(t)) 

Next, we evaluate the partial derivatives of r to 



n 



E,(t,E(-i)(t,E + 6))' 

t xx ^^(t.St-^t.S + ii)) 

" E x (t, £ + b)) E x (t, E(-D(t, £ + b)f 

-E,(t,S(- 1 )(t,E + 5)) + E t 

z x (t,z(-v(t,£ + b)) ' 
i 



Putting this together with (3.22), we end up with 

6 = --~a 2 t xx + I^ ^^EC- 1 )^ £ + &)) -, - _ ~ 

2 2 E x (t,EM)(t,E + 6)) 2 ffl H ' 1 ' " * 

Finally, recalling that H x = 1/a and E^ = l/o~, we can write the latter equation 

as 

b = l -a{t, E^t, £ + b)f E M (t, E ( - 1} (t, £ + 6)) + £ t (t, £ ( - 1} (i, £ + 6)) 

-^2^ 2 ^xx + £* 
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Moreover, since E 



and E, 



£f, this equation simplifies further to 



b = - l -a x {t, E^t, E + b)) + E t (t, &-V{t, t + b)) + l -a x - t t . 
The last equation holds due to the assumption that M and M are consistent. 



Thus, we have constructed a solution of the system (3.1), (3.2). 



Conversely, suppose that X, v, v form a solution of the system ( |3.1[ ), (3.2). 
Then, by the uniqueness result in Theorem [TJ the pair (X, v) has to coincide 
with the solution of (3.1) constructed in the proof of Theorem [TJ and the pair 
(X, v) has to coincide with the solution of ( |3.2[ ) constructed in the proof of 
Theorem [TJ In particular, it must hold 

(3.23) = J^— = -L. 

a(t,u) a{t,u) y/ a x 

One can now proceed as in the first part of the proof to deduce the existence of a 
function b : [0, T] — > R such that h is given by (3.16). Plugging this expression 



for h into (3.19) and proceeding as before, one shows that b must solve the 
ODE in Definition [2~J Mo reover, setting T = v o one easily verifies (3.4) 

and (3.5) by using (3.15) and E[t>(Jf (T))] = m , respectively. This shows the 
consistence of M and M. 



Finally, since the choice of v in the construction above was arbitrary among 
all functions in Q, we conclude that, if M and M are consistent, the system 
(3.1), (3.2) has uncountably infinitely many solutions. □ 



Remark 4. A careful reading of the proof of Proposition [7] shows that the system 
(3.1), (3.2) is solvable by a Brownian martingale X and functions v,v G Q 



under the weaker assumption that M satisfies the conditions of Theorem [TJ the 
martingale problem for | a(t, x) 2 4^ is well-posed and the diffusion coefficients 
a, a are continuously differentiate, bounded and bounded away from 0. The 
same is true for the statements in the upcoming Corollary [8j 

As a consequence of Proposition [7j we obtain the existence of couplings of 
consistent Brownian martingales. 

Corollary 8. Suppose that the Brownian martingales M and M satisfy the 
conditions of Theorem^ their diffusion coefficients are continuously differen- 
tiable and bounded away from ; and that M , M are consistent in the sense of 
Definition [2J Then, M and M can be defined on the same probability space to 
form a weak solution the degenerate system of SDEs 

(3.24) dM(t) = a(t, M{t)) dB(t), M(0) = m 

(3.25) dM{t) = a(t, M{t)) dB(t), M(0) = m 

on [0,T]. //, in addition, the function T in Definition^is such that x > T(x) 
for all x in the state space of M , then the inequality M{t) < M{t) holds for all 
t G [0,T] with probability 1. 
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Proof. Applying Proposition [7j we see that there is a probability space, a 
Brownian martingale X defined on this space and functions v,v G Q such 
that the equations j[3.1[ ), (3.2) hold. As we have seen in the course of the 
proof of Proposition [?Ttne processes E[t;(X(T))|J rX (-)], E[w(X(T))|J rX (-)] solve 



the degenerate system of SDEs (3.24), (3.25). Moreover, under the additional 



assumption of x > T(x) for all x in the state space of M, we have 

(3.26) E[v(X(T))\F x (.)\ = E[r(v(X(T)))\F x (.)\ < E[v(X(T))\F x (■)] 

with probability 1, which is the desired monotone coupling. 



□ 



Remark 5. In the setting of Remark [T] with a = 1, Remark [4] and Corollary [8] 
show that the degenerate system of SDEs 



(3.27) 
(3.28) 



dM(t) 
dM(t) 



dB(t), 
1 



q' T _ t (M(t)) 



dB(t) 



has a weak solution on [0, T] for any initial values mo, mo. We note hereby that 
the diffusion coefficient of M is continuous, bounded and bounded away from 
0, so that the martingale problem satisfied by M is well-posed. Moreover, since 
the weak solution constructed in the proof of Corollary [8] has the property that, 
for any t G [0, T], M(t) can be written as a deterministic function of M(t), 
the resulting process M is adapted to the filtration generated by the Brownian 



motion B. This shows that the SDE (3.28) has a strong solution for any initial 
value fho. We note that if the function 



Qt- 



[x 



(e x (t,q T -t(x))) x 



&xx{t,qT-t(x)) 

Q x (t,q T -t(x)) 



(loge x ) x (t,q T -t(x)) 



fails to be bounded, the existence of a strong solution does not follow from clas- 
sical existence theorems such as Theorem 2.9 in chapter 5 of [12]. In addition, 
in this case T = qo — D, where D is uniquely determined by 

rho = E[e(T,mo + B(T)+D)]. 

At this point, Corollary 8] shows that, whenever qo(x) < x + D holds for all 
x G M, the system (3.27), (3.28) has a weak solution on [0, T] starting from 
(mo, mo) such that the inequality m + B(t) < M(t) holds for all t G [0, T] with 
probability 1. 

To demonstrate the possible range of applications of Corollary [8j we give one 
immediate corollary. 

Corollary 9. Suppose that the Brownian martingales M and M satisfy the 
conditions of Theorem^ their diffusion coefficients are continuously differen- 
tiable and bounded away from 0, and that M, M are consistent in the sense 
of Definition [2j Assume further that the inequality x > T(x) holds for all x 
in the state space of M . Moreover, for any fixed < e < m , let r e , f e be the 
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first hitting times of the set [0, e] before T by the respective Brownian martin- 
gales (which we set to be equal to T if the set is not hit before T). Then, t € is 
stochastically dominated by f e in the sense that it holds 

(3.29) P(r e > t) < P(f £ > t) 

for all t G [0, T] . 

3.2. Systems of any number of equations. We now generalize the con- 
structions of the previous subsection to give a proof of Theorem [2] for a general 
index set A. 

Proof of Theorem [2J We assume first that the family M A , A G A is consistent 
in the sense of Definition [l] and will construct a solution of the system (1.2) 
of the desired type. Arguing as in the proof of Proposition [7j we see that 
it suffices to find a Brownian martingale X with a time-dependent generator 



1 „x 



(t,x] 



j^, as well as functions v x , 



A G A in Q such that the system of PDEs 



(3.30) 



u 



\ + a u xx 



u x V a 



x 



0. 



A G A 



A 



A G A 
A G A 

A G A. To this end, we fix a 



(3.31) 

(3.32) u X {T,-) = v . 

has a classical solution and K[v x (X(T))] = m x , 
A* G A and a function v x * G Q, and define X as the solution of (1.1) for the 
pair (M A *,t> A *), constructed in the proof of Theorem [TJ Then, the equations 
fl3~30| ), ( |33Tj ), ( |3~32| ) will hold for A*. Now, we set v x = (r^^)^ 1 ) o v x * for all 
A 7^ A* in A. Following the lines of the proof of Proposition [7], one checks that 
with this choice the equations (3.30), (3.31), (3.32) are satisfied for all A G A 
due to the consistency of the family M A , A G A. In addition, since there are 
uncountably infinitely many choices for the function v A * , there are uncountably 



infinitely many solutions of the system (1.2). 



Conversely, suppose that a Brownian martingale X and a family of functions 



A G A in Q solve the system (1.2). Then, by Proposition [7| every pair M Al , 
M A2 must be consistent in the sense of Definition |2j Therefore, the family M A , 
A G A is consistent in the sense of Definition [T] 

At this point, to show the statements (a) and (b) in the theorem, one only 
needs to follow the lines of the proof of Corollary [8] To prove statement (c), we 
apply Theorem 4.2 in [T7j to deduce that, under the assumptions in statement 
(c) in the theorem, there exists a version of the unique strong solution of the 
system (1.10), which is continuously differentiable in A for every t G [0, T], and, 
for all A G A, the derivative is a strong solution of the equation 



(3.33) 



dM> 



°a(A, £, M x (t)) + <Xj.(A, t, M x {t)) 



dB(t). 



Rewriting the latter equation in distributional form, we arrive at (1.11). □ 
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Remark 6 . A careful reading of the proof of Theorem [2] shows that the system 



(1.2) has a solution of the desired type and the statements (a), (b), (c) of 
Theorem [2] hold under the weaker assumption that there is a A* G A such 
that M x * satisfies the conditions of Theorem [TJ the martingale problems for 
\ a x (t, x) 2 A G A\{A*} are well-posed and the diffusion coefficients a x , 
A G A are continuously differentiable, bounded and bounded away from 0. 

Remark 7. Let A C R be an interval and define the Brownian martingales 
M A , A G A as in the second paragraph of Remark [TJ but choosing terminal 
conditions A (T, ■) depending on A. Letting q^_ t be the corresponding quantile 
functions, we recall from Remark [l] that the family M A , A G A of Brownian 
martingales with diffusion coefficients a x (t,x) := ^ A 1 y( x ) is consistent for any 



choice of initial values m A , A G A. In view of Remark [6] and Theorem [2] (a), the 
degenerate system of SDEs 

(3.34) dM(t) = dB(t), 

(3 - 35) iMHt) = WMMW) dB(tl AeA 

has a weak solution on [0, T] for any initial values mo and m A , A G A. We recall 
hereby from Remark [5] that, for every A G A, the martingale problem satisfied 
by M A is well-posed. Subsequently, we deduce as in Remark [5] that all processes 
M A , A G A are adapted to the filtration generated by the Brownian motion B 



by construction. This shows that the system (3.35) has a strong solution for 
any initial values m A , A G A. 

Next, define the constants D x , A G A by 



E[6 A (T, m + B(T) + D : 



If the functions T A := — D x decrease pointwise in A, then the functions 
pAi,A 2 ._ ^pA^(-i) Q r A 2 sat j s fy r Al,A2 (x) < x for all x in the state space of 
M Xi , whenever Ai < A2. In this case, we can conclude from Theorem^ (b) 



that the system (3.35) has a weak solution, for which M Al (i) < M 2 {t) holds 



for all t G [0, T] with probability 1, whenever Ai < A2. Finally, if the functions 
A 1 — y m x and (A, t, x) i-> ^ x satisfy the additional regularity assumptions 



of Theorem [2j (c), then there is a version of the growth process H defined there, 
which satisfies 

(3 - 36) - { wtm - wm^ Hx 1 ^ 

(3.37) ^°' A ) = ^ 

with £ being the distribution- valued Gaussian field in part (c) of Theorem [2] 
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Example 1. Based on Remark [7| we now give a numeric example of a situation, 
in which Theorem [2] (c) applies. To this end, we would like to choose the 
terminal conditions G A (T, •) as the functions 

(3.38) f x (x) := x 1{ X < } + X x 1{ X>0 } 

with A varying in A := [1, oo). However, the functions f x do not satisfy the 
differentiability assumption on G A (T, •) of Remark [Tj To avoid this problem, 
we choose a small smoothing parameter k G (0, 1) and set 

(3.39) 6 A (T,-) :=/ A *0 K , 

where <fi K is the normal density with mean and variance k. A straightforward 
computation then gives 

Q\t, x) = (0 A (T, ■) * <f> T _ t ) (x) = {f x * (f> T -t+ K ) (x) 

~ 27t exp ( - x 2 /(2(T - t + k))) +x + {X-1)x 9 T -t+ K {x), 

where (pT-t+K is the normal density with mean and variance T — t + k and 
^T-t+n is the corresponding cumulative distribution function. Therefore, the 
corresponding quantile functions q x , A 6 [1, oo) are given as the solutions of 



(3.40) x = (X- Df-^ e-^'^-^ + g x - t (x) 

+ (X-l)q X _ t (x) ^ t+K {q X _ t {x)). 
Differentiating both sides of this equation with respect to x, solving for . A 1 

\lT-t) 

and simplifying, we end up with the identity 

1 _ x - (A - l)y/(T - t + K)/(27r)e- g T-^) 2 /WT-t+-)) 

(3 ' 41) WW) " 4Jx) ' 

where the functions q x , X G [1,cxd) are given as the solutions of (3.40). We 
conclude from Remark [7] that the Brownian martingales M x , X G [1, oo) corre- 
spond!: 

values 



sponding to the diffusion coefficients cr x (t, x) = ^ A 1 y^ x y A G [1, oo) and initial 



m x = E[Q x (T,B(T))}= [ [ f x (y)4 >K (x-y)dyMx)dx 

Jr Jr 

f\y^ T+K (y)dy=(X-l)^=^, AG [1,00) 

V Z7T 

can be defined on the same probability space in such a way that the inequalities 
M Xl (t) < M X2 (t) are satisfied for all t G [0,T] and A x < A 2 in [l,oo) with 
probability 1. In Figure [TJ we demonstrate this fact by plotting sample paths 
of the Brownian martingales M 1 , M 2 , M 3 , M 4 ', M 5 for k = 10~ 4 obtained by 
using the same sample path of the driving Brownian motion. 
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4. Additional examples 

In this last section we demonstrate the results of Theorem [T] Proposition [6] 
and Theorem [2] on the example of the Kimura martingale 

(4.1) dM(t) = M(t)(l — M{t)) dB(t), t G [0,T], M(0) = m G (0, 1) 

and a time-changed version of it defined below. We recall that the state space 
of M is the open interval (0, 1) and that M does not exit from (0, 1) in finite 
time with probability 1 (see e.g. section 6 in |llj). 

4.1. Single equation. First, we fix a function v G Q, and seek a Brownian 
martingale X solving with M given by (4.1). To this end, we recall first 



that the transition densities of the Kimura martingale M can be written down 
explicitly (see [TT], [13]): 

a/2 



^(t,x;T,,) = (2.(T-t))-V 2 ||i_^ /2 



x exp — — r log 



2(T-t)V °x(l-y), 

i G [0,T), x,y E (0,1). In particular, the partial derivatives pf 1 , p^f , 
and exist and are continuous. In addition, a lengthy but straightforward 
calculation shows that 

(4.2) limp M (t, x; T, y) = limp M (t, ac; T,y) = 

and that the same is true for the partial derivatives pf 4 , p^f , and 
Therefore, we can apply Theorem [T] to conclude that there exists a unique 



Brownian martingale X, which satisfies (1.1) for these M and v. Moreover, 
Theorem [T] and its proof imply that X solves the SDE 

(4.3) dX(t) = h x (t, h^it, X(t))) /i ( - 1} (t, X(t))(l - h^it, X(t))) dB(t) 
on the time interval [0, T], where 

(4.4) h(t,x):= [ p M (t,x;T,y)v^ 1 \y)dy 

J {0,1) 



solves the problem (1.6) with a(t,x) = x 2 (l — x) 2 in the classical sense. 

To demonstrate an application of Proposition [6j we consider the process M 
obtained from the Kimura martingale M by the deterministic time change 

(4.5) a : [0, oo) -> [0,T), rnTfl-e" 8 ). 

In other words, M is the Brownian martingale given by the solution of the SDE 

(4.6) dM(t) = ^^£^dB(t),te[0,T), M(0) = m G(0,l). 

Now, by Proposition 5.22 (d) in chapter 5 of [T2], the almost sure limit M(T) : = 
lim^T M(t) = lim^oo M(s) exists and is given by 1 on a set of probability mo 
and by on a set of probability 1 — mo- Letting v = lr c ,oo) f° r some c G R, 
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we see from Proposition [6] that the equation (1.1) with M and this choice of 
v has uncountably infinitely many solutions, each of which satisfies an SDE of 
the form 

(4.7) dX W= (X W-c.Hc^ffl) dm 

(ci - coJVT - * 
for some constants Co < c < C\. 

4.2. Systems of equations. In order to give an example of an application of 
Theorem |2j we start with the Kimura martingale M and proceed as in Remark 
[T] to construct a family of Brownian martingales consistent with M . In this 
particular case, one has 

(4.8) S(t,x) = log " (1 - mo) , tf-V(t,x) 



[1 — x)m l — m Q + m e x 

(4.9) £ t (t, x) = 0, a x {t, x) = 1 - 2x. 

Hence, the backward heat equation of Remark [l] reads 

~ 1 ~ 1 mn e x+D + m n — 1 ~ 

4.io o = e t + - e xx + - -°— °— e x . 

2 2 m e x+u — m + 1 

Setting D = and letting 6(T, •) be a function in Q, we see that, for any fixed 
t G [0,T], it holds 0(t, ■) = 6(T, ■) * ipr-t with ^T-t being the appropriate 
transition density of the diffusion 

(4.11) dR(t) = ^°C +m °- 1 dt + dB(t). 

K J w 2 m e fl (*) - mo + 1 W 

For each t e [0, T], we write g^_ 4 for the spatial inverse of Q(t, •) and conclude 
from Remark [T] that the Kimura martingale M is consistent with the Brownian 
martingale solving the SDE 
(4.12) 

dM(t) = — L— dB(f), TO = E-re(r,iog 7 M(T)(1 - mo) v 



M(T))(to )/ 



Letting A C K be an interval, letting 0(T, ■) depend on A and proceeding 
as before, we obtain a family of consistent Brownian martingales M A , A G A. 
Moreover, the family M, M A , A G A satisfies the conditions described in Remark 
[6} Indeed, for every A G A, the boundedness, the boundedness away from 
and the continuity of the diffusion coefficients of M x shows that the martingale 
problem satisfied by M x is well-posed. Therefore, we may apply Theorem [2] (a) 
in this situation to conclude that the degenerate system of SDEs 

dM(t) = M(t)(l - M(t)) dB(t), 

dM A (t) = — — i— — di?(t), m x = E m ° [§ A ( T, log M{T) }\ ~ mo) 
1 ; {q x _ t )'{M\t)) yh L V 6 (1 - M(T))to 

has a weak solution on [0, T] . 
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If, in addition, the functions A (T, •) are increasing pointwise in A, then 
we can apply Theorem [2] (b) to conclude that there is a weak solution of the 
latter system of SDEs on [0,T] such that, for all Ai < A2 in A, the inequality 
M Xl (t) < M X2 (t) holds for all t G [0, T] with probability 1. 



If, in addition to the above, the function A 1— > m x is continuously differentiable 
and the function <r(A,t, x) := ^ x yj^ has continuous and bounded partial 



derivatives a\, a\\, a\ x , a x and a xx , then by Theorem|2j (c) there is a version of 
the growth process H(X,t) = M x (t), A G A, t G [0,T], which evolves according 
to 

(4,3) ^ A= f-[%V)]_ (^rm 

with £ being the distribution- valued Gaussian field in Theorem [2l part (c). 
Note that these equations coincide with the equations in Remark [7] with the 
difference being that here the functions q^-t are given by spatial inverses of 



convolutions with the heat kernel for the diffusion in (4.11), as opposed to the 
heat kernel of a standard Brownian motion in Remark!?! 
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